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淹没 水 射流 流 场 的 DES 模拟 与 实验 研究 


张 伟 康 灿 潘 寡 朱 洋 
(江苏 大 学 能 源 与 动力 工程 学 院 ， 镇 江 212013) 


摘 要 为 了 深入 分 析 淹 没 水 射流 流动 特征 ,综合 采用 数值 模拟 和 流 场 测量 进行 研究 ,数值 模拟 应 用 基于 S-A 模型 的 
DES 方法 , 并 通过 圆柱 绕 流 验 证 数值 方法 的 有 效 性 ; 流 场 测量 采用 粒子 图 像 测速 技术 . 研究 表明 : 采用 的 数值 模拟 方案 可 
以 准确 预测 圆柱 的 升力 系数 、 阻 力 系 数 和 尾 流 中 旋涡 脱落 的 斯 特 劳 哈 数 ; 模拟 和 实验 获得 的 射流 速度 沿 轴 向 和 径 向 的 衰减 
规律 吻合 ; 模拟 获得 的 射流 半 值 宽度 线 的 斜率 小 于 实验 结果 ,主要 是 由 于 对 剪 切 层 厚度 的 预测 不 足 所 致 ;， 在 垂直 于 射流 束 
模 截 面 上 , 涡 量 呈 现 双 峰 分 布 ; 射流 的 涡 量 随 着 射流 发 展 呈 总 体 衰减 趋势 . 
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DES Simulation and Experimental Investigation of Submerged Waterjet Flow 


ZHANG Wei KANG Can PAN Chen ZHU Yang 
(School of Energy and Power Engineering, Jiangsu University, Zhenjiang 212013, China) 


Abstract Numerical simulation and flow measurement are performed to explore the flow charac- 
teristics of submerged waterjet. The numerical simulation is based on the DES method and S-A 
model, and the validity of the numerical scheme is substantiated with flows around a circular cylinder. 
Particle image velocimetry is utilized to obtain flow velocity distribution. The results show that the 
lift and drag coefficients of the cylinder and the Strouhal number of the vortex shedding in the wake 
are accurately predicted with the numerical technique. Regarding the waterjet, numerically obtained 
jet velocity attenuation in streamwise and lateral directions agrees well with the experimental results. 
The approximate slope of the semi-width curve numerically obtained is smaller than its counterpart 
associated with experiment, as is mainly ascribed to the underestimation of the shear layer thickness 
with numerical approach. Dual-peak distribution of vorticity arises on the cross section of the waterjet 


stream. Vorticity magnitude decreases with the jet progression. 


Key words flow around a cylinder; submerged waterjet; detached eddy simulation; reynolds num- 
ber; nozzle 
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0 引 


淹没 水 射流 是 射流 的 重要 类 型 之 一 . 近年 来 , 淹 
没 水 射流 的 工程 应 用 范围 不 断 延 伸 。 从 目前 的 相关 
研究 来 看 , 淹没 水 射流 的 实验 研究 难度 较 大 , 一 是 采 
用 光学 手段 进行 流 场 测量 的 难度 较 大 ， 在 介质 中 撒 
播 的 示 足 粒子 对 介质 的 跟随 性 差 且 分 布 不 均匀 ; 二 


是 流 场 中 可 能 会 出 现 空 化 现象 ,对 于 空 化 发 生机 理 
的 认识 、 空 化 的 瞬 态 发 展 过 程 和 空 化 形态 的 捕捉 又 


属于 为 一 类 难题 。 从 淹没 水 射流 流 场 中 的 流动 特征 
来 看 ， 其 最 为 明显 的 特征 就 是 主流 与 静止 的 环境 水 
体 之 间 的 相互 作用 ， 苦 止 水 体 阻碍 主流 的 发 展 ， 而 
主流 与 环境 水 体 之 间 的 强 剪 切 作 用 又 可 能 入 生出 新 
的 流动 结构 . 从 完整 地 描述 淹没 水 射流 的 角度 出 发 ， 
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), 硕士 研究 生 , 主要 从 事 淹没 水 射流 流动 特性 和 气泡 运动 的 测量 方法 研究 . 


数值 模拟 能 够 发 挥 有 效 的 作用 . 

目前 的 灌流 数值 模拟 方法 主要 有 四 种 : 基于 雷 
诺 平 均 Navier-Stokes 方程 (RANS) 的 模拟 ， 直 接 
数值 模拟 (DNS)， 大 涡 模 拟 (LES) 以 及 分 离 涡 模 
拟 (DES)。 限于 计算 机 硬件 水 平 , DNS 方法 对 于 高 
雷诺 数 流动 的 数值 模拟 较 难 实现 . LES 方法 的 计算 
量 也 较 大 而 且 其 对 于 边界 层 发 展 和 分 离 的 模拟 欠 准 
确 . Spalart [1 在 1997 年 提出 了 DES 方法 , 其 基本 
思想 是 在 近 壁 区 采用 RANS 方法 , 用 注 流 模型 模拟 
其 中 的 小 尺度 脉动 运动 ， 在 远离 壁面 区 域 ， 将 济 流 
模型 耗 散 项 中 的 庙 流 尺 度 参数 与 一 常数 相 乘 ， 发 挥 
与 亚 格子 雷诺 应 力 模 型 相似 的 作用 ， 同 时 保证 了 近 
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壁 区 流动 和 远离 固体 壁面 的 区 域 的 模拟 允 近 物理 真 
实 。 从 对 大 尺度 流动 结构 的 捕捉 来 看 ，DES 方法 的 
适应 性 较 强 ,可 以 作为 对 实验 研究 的 有 效 补充 

本 文 首 先 设计 了 基于 Spalart-Allmaras 模型 的 
DES 数值 模拟 方案 , 在 雷诺 数 为 2x104 条 件 下 对 圆 
柱 绕 流 流 场 进行 模拟 以 验证 数值 模拟 方案 的 有 效 性 ， 
将 模拟 获得 的 升力 系数 、 阻 力 系 数 和 斯 托 劳 哈 尔 数 
与 实验 结果 进行 对 比分 析 . 进而 采用 DES 方法 对 淹 
没 水 射流 流 场 进行 模拟 ， 模 拟 中 综合 考虑 射流 的 速 
度 衰 减 、 半 值 宽 、 涡 量 分 布 等 表征 射流 特征 的 重要 
因素 ; 同时 采用 粒子 图 像 速 度 场 仪 辅 以 效 光 粒子 对 
淹没 水 射流 流 场 进行 测量 ,结合 数值 模拟 和 实验 结 
果 对 淹没 水 射流 流动 进行 系统 性 的 解释 ， 并 对 比 数 
值 模拟 与 实验 结果 之 间 的 异同 . 


1 数值 方法 
1.1 流体 控制 方程 


对 于 带 有 固体 边界 的 非 定常 理性 流动 ， 采 用 
Navier-Stokes 方程 作为 控制 方程 , 其 积分 形式 为 : 


D 了 二 

一 Wd Fe 一 上 vjd9=0 1 

ot | 4 ( ) 外 
式 中 , W, 天 ,及 分 别 为 守恒 矢量 、 对 流通 量 和 黏 性 
通 量 , 其 具体 含义 见 文献 [2]， 
1.2 Spalart-Allmaras 模型 

S-A 模型 以 涡 黏 系数 为 变量 ， 它 的 积分 形式 

为 : 


了 于/ oan+ (Wer Bras= QrdD (2) 
Ot jp Ele) 0 


i 

QT=(cp1S + A Vcufs (3) 十 (A0) (5) 
分 别 为 汕 流 涡 黏 系数 的 对 流通 量 、 黏 性 通 量 和 源 项 ， 
其 具体 含义 以 及 计算 常数 参见 文献 [3]。 
1.3 DES 满 流 模型 

Spalart 通过 修改 Spalart-Allmaras 模型 中 的 长 
度量 4 提出 了 基于 S-A 模型 的 DES 方法 . 其 微分 形 
式 为 : 


Ds _ ~ 了]2 
元 = cp15% 一 cw1 fw 四 站 


=[V . ((v + td)VD) + coz (VD 
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式 中 ，d 为 空间 点 到 固体 壁面 的 最 短 距 离 ，5 为 
修改 后 的 涡 量 表达 式 . DES 方法 将 上 式 中 的 d 用 
d=(d, CppsA) 来 代替 .其 中 A 为 当地 网 格 的 最 大 
尺寸 , Cbas 为 常数 , 取 0.65. 显然 , 在 固体 壁面 附近 
有 d<Cpgs, 于 是 d= d, 这 时 候 的 DES 模型 就 是 式 
(4) 所 描述 的 S-A 模型 . 当 d>CpEs 时 , d = CppsAA， 
此 时 清流 涡 黏 系数 的 训 减 由 当地 的 网 格 尺度 来 确定 . 
通过 式 (4) 可 以 得 知 , 当 汕 流 涡 务 系数 的 源 项 和 汕 流 
涡 黏 系数 的 衰减 项 达到 平衡 时 , 5 正比 于 5d2, 如 果 
d= CpgsA, 则 1 正比 于 $A?, 此 即 亚 格 子 模型 。 

网 格 布置 是 影响 DES 求解 结果 准确 性 的 关键 . 
近 壁 区 网 格 的 最 大 间距 和 网 格 点 到 壁面 的 距离 决定 
了 涡 黏 性 模型 与 亚 格子 模型 的 交界 面 位 置 ， 如 果 璧 
面 附近 过 度 加 密 ， 就 会 过 早 地 引入 亚 格子 模型 ， 因 
此 对 边界 层 的 分 辨 尤为 重要 .通过 引入 参数 ra 确 
保 边 界 层 内 为 涡 黏 性 模型 ,将 长 度 尺 寸 4 关系 式 改 
写 为 


d= qd— famax(0,d— CprsA) (7) 
其 中 : 
万 =1 一 tanh((8ra)3) (8) 
WV (9) 
a A Ui jUi,ir?2d 


fa 在 边界 层 内 为 0, 在 其 他 区 域 等 于 1, vi 为 运动 涡 
黏 性, v 为 运动 分 子 黏 性 , Ui,; 为 速度 梯度 , r=0.41， 
通过 函数 关系 实现 亚 格子 应 力 模 型 和 涡 务 性 模型 的 
切换 . 


2 数值 算 例 验证 


圆柱 绕 流 是 验证 数值 模拟 的 典型 流动 。 李 燕 玲 
等 向 以 及 尹 纪 富 等 回 都 采用 DES 方法 对 圆柱 绕 
流 问 题 进行 了 模拟 ， 并 且 都 取得 了 与 实验 数据 相近 
的 结果 . 此 处 为 校 验 本 文 所 设计 的 DES 方案 的 有 效 
性 ,选取 雷诺 数 为 2x104 的 圆柱 绕 流 为 对 象 ， 模 拟 
圆柱 绕 流 时 出 现 的 非 定常 流动 特征 。 
2.1 圆柱 绕 流 验证 

1) 网 格 划 分 和 边界 条 件 

坐标 系 定 义 如 图 1 所 示 : 坐标 原点 为 圆柱 轴线 
的 中 点 ，X 为 来 流速 度 方 向 ，2 为 圆柱 轴线 方向 
以 圆柱 直径 D 为 特征 尺度 ， 整 个 长 方 体 计算 域 的 
尺寸 为 : -15D < XX 过 35D, -10D < YY <10D 和 
一 5D < 2 < 5D, 如 图 1 所 示 。 该 计算 域 可 以 保证 来 
进行 网 格 加 密 , 以 圆柱 的 轴 截 面 为 中 心 , 在 8D x 8D 
范围 内 沿 径 向 按 比例 逐步 放大 网 格 , 如 图 2 所 示 . 计 
算 域 的 网 格 总 数 约 为 159 万 。 
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所 有 的 固体 壁面 处 采用 无 滑 移 边界 条 件 ， 计 算 
域 的 入 口 设 定 均 匀 来 流 条 件 ， 出 口 采用 压力 出 口 边 
界 条 件 , 出 口 面 上 的 压力 梯度 设 为 零 . 以 均匀 进口 速 
度 分 布 作 为 非 定常 计算 的 初始 条 件 ， 采 用 SIMPLE 
压力 速度 耦合 欠 代 方法 , 压力 项 离散 采用 二 阶 精 度 ， 
根据 雷诺 数 的 量 级 , 将 时 间 步 长 设 为 0.025 s。 


图 1 圆柱 绕 流 计算 域 模型 
Fig. 1 Computational domain for flows around a circular 
cylinder 


Fig. 2 Mesh deployment on the cross section vertical to 
cylinder axis 


2) 圆柱 绕 流 计算 结果 分 析 

阻力 系数 和 升力 系数 反映 了 流体 作用 在 圆柱 表 
面 上 的 载荷 ， 而 斯 特 劳 哈 数 则 用 来 描述 圆柱 尾 流 中 
旋涡 的 脱落 ,该 3 个 无 量 纲 数 对 于 圆柱 绕 流 问题 最 
具 代表 性 . 在 表 1 中 , 将 数值 计算 得 到 的 时 均 阻力 
系数 Ca、 升力 系数 Cl 和 斯 特 劳 哈 数 54 与 文献 中 的 
实验 结果 和 模拟 结果 进行 了 对 比 .。 可 以 看 出 ， 和 其 
他 方法 相 比 , 本 文 的 模拟 结果 更 接近 实验 结果 。 

尾 流 中 旋涡 的 周期 性 脱落 是 圆柱 绕 流 的 特征 现 
象 ， 随 着 旋涡 的 脱落 ， 圆 柱 所 受到 的 阻力 和 升力 会 
出 现 周 期 性 变化 . 图 3 中 表示 出 了 垂直 于 圆柱 轴 截 
面 上 的 瞬时 流 场 结构 ,分 别 以 压力 等 值 线 和 速度 分 
布 表示 . 可 以 看 出 , 尾 流 中 形成 两 列 上 、 下 交替 脱落 
的 旋涡 ,离开 圆柱 的 流向 距离 越 大 ,旋涡 强度 越 弱 ; 


从 速度 分 布 的 形态 来 看 , 尾 流 基本 上 呈现 振荡 状态 ， 
且 越 靠近 圆柱 ,外 部 流动 与 尾 流 间 的 速度 梯度 越 大 ， 
这 也 解释 了 近 圆 柱 区 下 游 襄 基 值 高 的 原因 . 


表 1 圆柱 绕 流 参数 对 比 
Table 1 Comparison of parameters associated 


with flows around a circular cylinder 


Re=2x104 Cd 时 均值 C1 幅 值 St 
Lu [6] (RANS) 1.15 0.60 0.20 
Sun(DNS) 1.30 0.80 0.29 
Yokuda [ (实验 ) 1.20 一 0.21 
苏 铭 德 [8] (LES) 1.22 0.80 0.23 
本 文 结 果 1.18 0.62 0.22 


(a) 压力 分 布 


(a) Pressure distribution 


(b) 速度 分 布 
(b) Velocity Instantaneous distribution 
图 3 瞬时 流 场 分 布 


Fig. 3 Instantaneous flow characteristics 


为 深入 理解 圆柱 绕 流 中 涡 结构 的 三 维 形态 ， 提 
取 了 瞬时 涡 量 等 值 面 ,并 从 图 4 所 示 的 视角 进行 观 
察 。 涡 量 等 值 面 表达 了 涡 结构 之 间 的 慕 连 、 合 并 与 
分 离 ， 并 有 反映 出 涡 结构 的 尺度 与 涡 强度 。 从 图 4 可 
以 看 出 , 在 主流 方向 上 存在 着 交替 发 展 的 涡 街 现象 ， 
沿 圆柱 轴 向 还 出 现 涡 结构 间 的 相互 交叉 .该 三 维特 
征 与 已 发 表 的 实验 和 模拟 结果 相近 1。 随 着 主流 
向 下 游 发 展 ， 涡 结构 的 整体 强度 减弱 ， 尤 其 是 上 下 
两 端的 涡 结构 逐渐 崩塌 ， 圆 柱 绕 流 的 二 维 假设 不 复 
成 立 。 从 对 大 不 度 流动 结构 的 捕捉 和 对 流 场 宏观 参 
数 的 求解 来 看 , 本文 的 DES 方案 具有 明显 的 优势 . 
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图 4 三 维 瞬时 涡 量 图 


Fig. 4 Three-dimensional instantaneous distribution pattern of iso-vorticity elements 


2.2 淹没 水 射流 的 数值 模拟 方案 

采用 经 过 验证 的 DES 方案 对 淹没 水 射流 进行 
数值 模拟 . 采用 的 射流 喷嘴 出 口 直径 为 D=2.0 mm， 
整个 计算 域 设 为 长 方 体 , 便于 直角 坐标 系 的 应 用 ; 射 
流 介 质 为 常温 水 . 在 计算 域 入 口 设置 速度 边界 条 件 ， 
设 定 喷嘴 的 入 口 速度 分 别 为 7.0、9.5 和 18.0 m/s。 
计算 域 的 直角 坐标 系 定 义 如 图 5 所 示 : 将 喷嘴 出 口 
中 心 定义 为 坐标 原点 ，2 方 疝 为 主流 方向 ,，X、Y 
为 垂直 于 来 流 的 速度 方向 。 环境 水 体 的 计算 域 尺 
寸 大 小 为 : -15 mm<X<15 mm, 一 15 mm<Y<15 
mm, -15 mm<2Z<100 mm, 对 计算 域 网 格 进行 局 部 
加 密 , 如 图 6 所 示 . 计算 域 的 网 格 总 数 约 为 128 万 。 


30 100 


图 5 淹没 水 射流 计算 域 几何 模型 


Fig. 5 Computational domain for submerged waterjet 


图 6 计算 域 截面 网 格 


Fig. 6 Cross-sectional view of meshes discretizing the 


computational domain 


以 进口 边界 条 件 作 为 计算 的 初始 条 件 , 采用 SIMPLE 
压力 速度 确 合 迭代 方法 ,压力 项 离散 和 采用 二 阶 精 度 ， 
时 间 步 长 设 为 0.0001 s。 


3 淹没 水 射流 光学 测试 方案 


为 了 进一步 评估 DES 数值 模拟 结果 的 有 效 性 并 
补充 射流 信息 , 采用 粒子 图 像 测 速 技术 (PIV) 8 对 
淹没 水 射流 流 场 进行 了 实验 研究 。 根据 射流 压力 与 
喷嘴 出 口 流 速 之 间 的 关系 ， 在 实验 中 设 定 与 数值 模 
拟 相 同 的 三 种 射流 速度 : 7.0、9.5 与 18 m/s. 淹没 射 
流 实验 中 , 由 于 射流 与 周围 水 体 间 的 速度 梯度 较 大 ， 
可 能 会 产生 空 化 现象 ， 示 踪 粒 子 的 物性 和 播 洒 方案 
极为 关键 . 本 次 实验 中 选用 罗丹 明 B 作为 示 踪 粒子 ， 
罗丹 明 B 对 光线 最 佳 吸 收 和 放射 的 波长 分 别 为 526 
nm 和 580 nm， 能 合适 的 吸收 光谱 ， 可 以 很 好 的 配 
合 激光 设备 的 使 用 . 选取 的 罗丹 明 B 的 密度 为 0.79 
g/mL, 在 水 中 的 溶解 度 为 1 mg/mL, 配置 水 溶液 浓 
度 为 15 mg/mL. PIV 测试 系统 的 基本 原理 如 图 7 所 
示 。 与 一 般 实 验 不 同 的 是 ， 此 处 直接 在 被 测 流 场 中 
均匀 散布 示 踪 粒子 罗丹 明 B， 因 为 通过 实验 的 内 循 
环 系统 , 粒子 能 很 快 地 分 散 均 匀 , 当 流 场 稳定 后 , 利 
用 偏光 镜 组 将 脉冲 激光 转化 为 脉冲 片 光 ， 入 射 至 需 
要 测试 的 流 场 区 域 ,然后 通过 高 速 数码 相机 进行 双 
帧 双 曝 光 图 像 采 集 , 记录 该 时 刻 的 粒子 图 像 并 储存 ， 
最 后 通过 图 像 处理 软 件 ， 采 用 互相 关 算 法 进行 图 像 
处 理 , 获得 瞬时 流 场 的 速度 分 布 。 

实验 中 采用 Lavision 公司 生产 的 PIV 系统 , 示 
踪 粒 子 直 径 为 20~50 hm。 图 像 采 集 系统 主要 包括 
高 速 数码 相机 、 双 脉冲 激光 器 、 导 光 及 偏光 镜 组 部 
件 、 同 步 器 及 装 有 采集 系统 软件 的 计算 机 . CCD 相 
机 型 号 为 InagerProSX 5M， 其 分 辨 率 为 2456x2058 
pixel， 像 素 尺 二 为 3.45 hmx3.45 um， 最 大 帧 频 14.5 
Hz， 双 帧 双 上 曝光 模式 下 ， 最 小 帧 之 间 的 时 间 间 隔 为 
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600 ns， 镜头 采 用 Nikon 型 号 为 AF Nikkor 50 mm 
f/1.8 D 的 相机 镜头 ;激光 器 采用 Nd:YAG 双 脉 冲 
激光 器 ， 最 高 频率 15 Hz， 最 大 能 量 70 mJ, 光波 长 
为 532 nm。 图 像 采 集 及 处 理 采 用 DaVis 和 Tecplot 
软件 . 


双 脉 冲 激光 器 


计算 机 
图 7 PIV 测试 系统 示意 图 


Fig. 7 Schematic view of the PIV measurement system 
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现 径 向 速度 衰减 曲线 基本 重合 . 因此 ， 淹没 射流 径 
向 速度 衰减 规律 可 用 高 斯 分 布 描述 为 : 


U/Um = exp(—0.74(z/y12)") (10) 


o DES,W0=7 m/s 
07 | 0o DES,W=9.5 m/s 
A DES,W=18 m/s 

0.6 上 | 一 一 实验 ,WW=7 m/s 

| | - - -实验 ,WW=9.5 m/s 
05 上 |-…… 实 验 ,WW=18 m/s ， 
| | -.- 理 论 S 


TT | TT LI 
8 10 12 14 16 18 20 22 24 
z/D 


图 8 轴 向 速度 的 沿 程 衰减 


< Fig. 8 Attenuation of axial velocity in streamwise direction 
©O 

局 。 4 结果 分 析 与 讨论 

ee 1.0 0 DES,Z2=50 mm 
O ee | :DE Ze 
对 于 非 淹没 水 射流 ， 当 射流 进入 自 相似 区 后 ， 08 上 一 和 0 
元。 无量 纲 化 的 速度 剖面 与 雷诺 数 和 距 喷嘴 的 距离 无 .| se 
Ql 关 83。 而 对 于 淹没 水 射流 ， 尚 无 相关 的 报道 Raos| 

一 图 8 所 示 为 射流 轴线 上 的 无 量 纲 速度 分 布 . 在 全 04[ 

一 距 喷嘴 最 近 的 一 段 轴 向 距离 (2/D < 10) 内 , 由 于 喷 of 

p> 嘴 出 口 的 空 化 现象 和 流 束 表面 的 反光 现象 ，PIV 测 02 上 U/U, e040) 

= 速 时 的 采样 规律 性 差 , 无 法 获得 稳定 的 采样 信号 , 所 IE 

王 以 实验 结果 从 DD > 10 开始 . 从 模拟 结 采 来 看 ， 二 0 0.2 0.4 06 0.8 1.0 1.2 14 16 18 2.0 
© 条 轴 向 速度 衰减 曲线 基本 重合 ， 轴 向 速度 衰减 的 相 x/yin 


似 规律 明显 . 在 近 喷 嘴 区 , 模拟 结果 中 轴线 上 的 射流 
速度 大 小 基本 保持 不 变 , 该 段 为 射流 的 初始 段 , 射流 
的 集束 性 强 . 实验 获得 的 三 条 曲线 不 光滑 是 由 于 环 
境 水 体 对 射流 的 扰动 已 影响 到 射流 轴线 ， 这 一 点 与 
非 淹没 射流 有 着 明显 的 区 别 ia. 尽管 如 此 ,三 条 实 
验 曲 线 的 形态 并 无 明显 的 突变 ， 并且 与 淹没 射流 轴 
向 速度 衰减 理论 曲线 0 吻合 较 好 .从 模拟 结果 与 
实验 结果 的 对 比 来 看 ， 速 度 沿 流向 的 变化 趋势 较 吻 
合 , 证 明了 淹没 水 射流 的 能 量 的 沿 程 衰减 规律 . 

图 9 所 示 为 射流 初始 速度 18 m/s 时 ， 三 种 
实验 工 况 以 及 对 应 的 三 组 数值 模拟 在 距离 喷嘴 50 
mm、60 mm 和 80 mm 处 径 向 速度 的 无 量 纲 化 分 布 . 
可 以 发 现实 验 与 模拟 的 结果 高 度 吻 合 ， 射 流 各 断面 
上 轴 向 流速 的 径 向 分 布 具有 很 强 的 自 相 似 性 . 并 且 
通过 与 已 发 表 文献 [15] 中 的 计算 结果 进行 对 比 ,发 


图 9 径 向 流速 的 无 量 纲 化 分 布 


Fig. 9 Velocity distribution in lateral direction 


2) 淹没 水 射流 的 半 值 宽 特征 

图 10 所 示 为 淹没 水 射流 的 无 量 纲 半 值 宽度 ( 射 
流 轴 向 流速 最 大 值 的 一 半点 距离 轴线 的 径 向 距离 ) 分 
布 . 在 不 同 压力 条 件 下 , 模拟 获得 的 三 条 半 值 宽 线 基 
本 重合 , 且 在 近 喷 嘴 区 (2Z/D < 10) 半 值 宽大 小 基本 
保持 不 变 , 而 实验 获得 的 三 条 曲线 也 基本 重 乔 . 然而 
实验 和 数值 模拟 结果 之 间 存 在 一 定 的 偏差 。 半 值 宽 
斜率 的 大 小 表征 了 射流 径 向 的 扩展 程度 ,根据 前 人 
的 研究 结果 ， 淹 没 射流 半 值 宽 斜率 大 小 约 为 0.095. 
若 对 模拟 和 实验 的 结果 进行 线性 拟 合 ， 则 实验 数据 
线 的 斜率 明显 大 于 模拟 结果 ， 且 实验 的 速度 半 宽 扩 
展 率 近 似 为 0.095. 从 速度 分 布 的 角度 来 看 ， 射 流速 
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度 自 轴线 沿 径 向 逐渐 减 小 ， 在 模拟 结 采 中 反映 这 一 
趋势 较为 平缓 旦 速度 梯度 较 小 ， 在 实验 结果 中 , 在 
射流 束 与 环境 水 体 交 界面 附近 ， 出 现 剧烈 的 速度 下 
降 即 较 高 的 速度 梯度 ， 所 以 导致 图 10 所 示 的 差异 . 
从 男 一 角度 考虑 ， 数 值 模拟 中 对 于 淹没 水 射流 受到 
环境 水 体 的 阻力 预测 不 足 . 


—0O—DES,W=7 m/s 

> | —O—DES, W=9.5 m/s 
上 | DES, Wi=18 m/s 
一 本 实验 ,=7 m/s 
一 < 实验 ,=9.5 m/s 
一 -实验 ,1W=18 m/s 


0 5 10 14 20 25 30 35 40 
z/D 


图 10 半 值 宽 的 无 量 纲 化 处 理 


Fig. 10 Variation of semi-width with jet development 


3) 射流 场 中 的 涡 量 分 布 

图 11 所 示 为 射流 初始 速度 18 m/s 时 数值 模拟 
和 实验 获得 的 径 向 涡 量 分 布 . 图 中 以 喷嘴 直径 为 轴 
向 距离 的 无 量 纲 基准 ， 以 径 向 涡 量 最 大 值 为 纵 坐 标 
的 无 量 纲 化 基准 .从 涡 量 图 可 以 看 出 ， 在 不 同 截面 
上 DES 数值 模拟 与 实验 的 平均 涡 量 分 布 相近 , 在 射 
流 轴线 上 , 涡 量 的 绝对 值 近似 为 0, 自 射流 轴线 沿 径 
向 至 射流 两 人 出 边界 ， 涡 量 先 增 大 后 减 小 ， 并且 轴 比 
两 侧 的 涡 基 值 具 有 明显 的 对 称 性 . 这 是 因为 轴线 为 
射流 的 核心 区 , 轴 向 射流 主要 受到 水 流 的 阻力 影响 ， 
受到 的 两 侧 水 体 扰 动作 用 很 小 ， 而 在 核心 区 外 的 射 
流 与 周围 水 体 间 存 在 较 大 的 速度 梯度 ， 从 而 产生 较 
高 的 涡 量 ， 因 此 射流 与 环境 水 体 间 的 剪 切 和 摊 混 导 
致 剧烈 的 能 量 传递 与 耗 散 . 

图 12 所 示 为 PIV 实验 获得 的 瞬时 涡 量 分 布 , 在 
射流 上 下 两 侧 涡 量 矢量 方向 相反 ， 上 侧 为 顺 时 针 方 
向 , 下 侧 为 逆 时 针 方向 , 并且 其 大 小 呈 对 称 性 分 布 ， 
轴线 上 涡 基 值 近似 为 0. 在 近 喷 嘴 区 , 涡 量 值 最 大 且 
较 稳 定 ， 这 是 因为 在 射流 的 初始 段 ， 射 流 束 保持 了 
较 高 的 初始 动能 , 集束 性 强 . 随 着 射流 的 发 展 , 射流 
与 周围 水 体 发 生 相 互 作用 ， 能 量 逐 渐 耗 散 ， 同 时 油 
基 值 沿 流向 减 小 上 下 涡 对 之 间 的 时 序 性 减弱 ， 上 自 
涡 量 集中 区 的 边缘 入 生出 更 小 尺度 的 涡 ， 直 观 的 反 
映 出 了 淹没 水 射流 的 涡 量 变化 特征 。 从 滴 量 分 布 的 
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演化 过 程 也 可 以 推 知 ,淹没 水 射流 空 化 的 发 生 与 涡 
量 密切 相关 ， 只 可 能 出 现在 近 喷 嘴 区 . 


o DES,Z10mm 
oOo DES, Z-30 mm 
A DES,Z2=50mm 

一 一 实验, 2=10 mm 


=- - -实验 , Z-30 mm 
实验 , Z=50 mm 


/Om 


2 EE PE PE SE OY SE OE PS PE EE PE EE ES | 
-3.5-3.0-2.5-2.0-1.5-1.0-05 0.0 05 1.0 15 2.0 2.5 3.0 3.5 
x/D 

图 11 射流 场 中 的 涡 量 分 布 


Fig. 11 Cross-sectional vorticity distribution 


60 80 


40 
Z/mm 
图 12 瞬时 涡 量 分 布 


Fig. 12 ‘Temporal vorticity distribution 


5 结 论 


1) 通过 高 雷诺 数 圆柱 绕 流 验 证 了 DES 数值 方 
案 的 有 效 性 ， 所 预测 的 升力 系数 、 阻 力 系数 和 斯 特 
劳 哈 数 较 已 有 文献 结果 更 接近 于 实验 测量 值 ， 并且 
DES 数值 模拟 绪 得 了 圆柱 绕 流 的 展 向 流动 结构 形态 
及 非 均匀 特征 . 

2) DES 数值 模拟 得 到 的 淹没 水 射流 中 轴 向 速度 
和 径 向 速度 分 布 规律 与 PIV 实验 获得 的 结果 吻合 度 
高 . 对 于 淹没 水 射流 的 半 值 宽 分 布 规律 ,数值 模拟 
结果 较 实验 结 采 的 斜率 小 ， 说 明 射 流 与 环境 流体 间 
的 剪 切 效应 剧烈 ,数值 方法 难以 准确 预测 剪 切 层 厚 
度 和 速度 梯度 . 

3) 从 涡 量 分 布 来 看 ,射流 束 轴 心 线 两 侧 涡 量 值 
呈 轴 对 称 分 布 , 射流 束 轴 心 线 处 出 现 涡 量 的 最 小 值 . 
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随 着 射流 向 下 游 发 展 ， 流 束 不 断 沿 径 向 扩展 ， 涡 量 
总 体 呈 衰减 趋势 ， 空 化 发 生 的 可 能 性 降低 。 
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